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1«   Introduction. 

This  note  is  concerned  with  an  axisymmetric  disturbance 
which  travels  over  the  surface  of  a  sphere  and  obeys  the  equation 

n)  ^J:-..^  sin  Q  21  =  J_£i 

^^'  sin  Q  9£?  ^^"  ■=*  3©    ^2  3^2 

or,  if  cos  ©  =  X,  the  equation 

(2)        5l  (1-^)11  =  ^^  • 

The  disturbance  is  supposed  to  be  generated  in  the  neighborhood 
of  a  point  P  on  a  sphere;  and  the  object  is  to  trace  the  distur- 
bance as  it  moves  between  P  and  the  antipodal  point. 

The  equation  (1)  may  be  regarded  as  a  specialization  of  the 
wave  equation 

Where,  in  spherical  coordinates,  ©  is  the  co-latitude,  ({>  is  the 
longitude,  and  p  is  the  radial  distance,  while  t  is  the  time. 
There  are  nxomerous  problems  in  mathematical  physics  for  which  a 
solution  of  (3)  which  is  independent  of  p  and  (j)  has  interest. 
The  equation  (2)  arises  in  connection  with  the  analysis  of  sound 
waves  in  a  narrovr  region  bounded  by  tvxo  concentric  spheres, 
elastic  waves  in  a  spherical  tightly-stretched  membrane,  water 
waves  on  the  surface  of  a  shallow  spherical  ocean  and  many 
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other  problems  concerned  with  wave  motions  on  vibrations.   Con- 
sequently, the  equation,  subject  to  appropriate  initial  and 
boundary  conditions,  certainly  has  been  solved  more  times  than 
we  are  avrare  of.   However,  we  have  not  seen  an  analysis  of  the 
problem  under  consideration;  and  In  all  other  analyses  of  (2) 
that  we  have  seen  the  solution  is  presented  in  the  form  of  an 
infinite  series  of  eigenfunctions ,   This  kind  of  representation 
is  of  course  quite  satisfactory  for  some  purposes  but  it  throws 
little  light  on  the  intrinsic  nature  of  the  solution  when  we 
wish  to  consider  it  as  representing  a  progressive  disturbance. 
The  purpose  of  what  follows  Is  to  show  that  a  certain  disturb- 
ance which  satisfies  (1)  can  be  described  in  such  a  way  that  the 
description  is  similar  to  the  D'Alembert  description  of  the 
progress  of  a  disturbance  in  a  stretched  string. 

The  passage  of  a  x-jave  over  the  entire  surface  of  the  earth 
is  certainly  not  an  impossible  phenomenon.   It  is  probable  that 
in  prehistoric  times  there  were  some  terrific  submarine  volcanic 
explosions  which  created  water  waves  of  such  force  that  they 
traveled  enormous  distances;  and  it  is  even  possible  that  some 
of  these  waves  were  propagated  back  to  the  point  of  disturbance 
after  collision  at  the  antipodal  point.   This  is  just  the  be- 
havior which  certain  atmospheric  waves  exhibited  during  the 
eruption  of  Krakatoa,  the  most  gigantic  volcanic  action  known  to 
man, 

Krakatoa  is  the  anglicized  name  of  a  volcano  in  the  Sundra 
Strait,   Situated  between  Java  and  Sumatra;  and  lying  at  the 
intersection  of  submarine  volcanic  fissures,  Krakatoa  showed 
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signs  of  uneasy  activity  in  the  spring  and  summer  of  I883.   It 
exploded  on  Aug.  26,  I883;  and  violent  paroxysms,  reaching  a 
peak  on  the  27th,  continued  for  two  days.   Neighboring  Islands 
were  either  devastated  or  annihilated.   Approximately 
200,000,000,000  cubic  feet  of  material  were  extruded;  and  it  has 
been  estimated  that  the  energy  of  the  excited  atmospheric  wave 
alone  x-jas  some  one  thousand  times  that  of  the  energy  released  by 
the  atomic  bomb  explosion  at  Nagasaki.   The  subsequent  effects 
of  the  eruption  were  startling;  some  were  frightful,  some  were 
wonderful;  all  were  almost  incredibly  awesome,  and  many  were 
world  wide.   [1] " 

Of  all  the  numerous  spectacular  phenomena  excited  by  the 
eruption,  the  global  behavior  of  the  atmospheric  pressure  viave 
of  August  27,  1883  is  the  one  of  most  interest  in  the  present 
context.   This  wave  was  created  by  the  most  intense  sudden  ex- 
plosion of  the  eruption.  At  some  stations  at  least  seven 
passages  of  the  pressure  vjave  i-jere   recorded,  four  as  the  wave 
came  from  Krakatoa;  and  three  as  the  wave  came  from  the  anti- 
podal reflection  point  located  in  Venezuela.   For  each  of  the 
first  two  passages,  recording  barometers  at  many  stations  showed 
that  as  the  wave  arrived  and  passed  over,  the  pressure  was 
suddenly  increased  from  p  to  a  viavy  summit  then  sharply  decreased 
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I  See  [1]:-  Eruption  of  Krakatoa  and  Subsequent  Phenomena. 
Report  of  the  ICrakatoa  Committee  of  the  Royal  Society  of 
London.   (I888). 

The  author  wishes  to  thank  Prof,  Stoker  for  calling  his 
attention  to  Krakatoa,  and  for  the  helpful  conversations  which 
led  to  the  following  analysis. 
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to  a  lower  level  of  values  oscillating  about  p-k  (k  >  0) ;  and 
that  the  pressure  variation  vjas  detected  for  nearly  two  hours. 
In  subsequent  passages  the  wave  became  deformed  and  the  pressure 
variation  lost  the  character  described  above  although  the  return- 
ing wave  still  caused  rises  in  pressure. 

The  following  dlagraras,  taken  from  [1],  show  the  progress 
of  the  vrave  duj?ing  the  first  four  envelopments  of  the  earth. 
The  circles  on  the  right  represent  the  hemisphere  of  the  earth 
vjhich  has  Krakatoa  as  its  pole;  those  on  the  left  represent  the 
opposite  hemisphere  with  the  Antipodes  as  pole.   The  vertical 
line  through  k  is  the  meridian  which  passes  through  the  North 
pole.  The   asterisk  in  the  first  left  hand  circle  shows  the 
position  of  Cape  Horn  at  the  southernmost  tip  of  South  America. 
The  contour  lines  show  the  progress  of  a  phase  of  the  pressure 
vjave  at  each  successive  even  hour  of  GreenX'Jich  mean  time  begin- 
ning at  [{.  hrs.  27  Aug.  I883,   It  will  be  noticed  that  the  lines 
representing  the  successive  positions  depart  more  and  more  from 
circular  form  as  the  elapsed  time  from  the  initiation  of  the 
disturbance  increases.   This  indicates  that  the  velocity  of 
various  points  of  the  wave  front  v;as  not  constant. 


Wave  Ko.  1 
First  Passage 
From  Krakatoa  to  Antipodes 


Wave  No,  2 
First  Passage 
From  Antipodes  Back  to  Krakatoa 


I.    . 


Wave  No.  3 
Second  Passage 
Krakatoa  to  Antipodes 


Wave  No,  ii 

Second  Passage 

Antipodes  Back  to  Krakatoa 


2»   An  Interpretation 
In  order  to  have  a  concrete  interpretation  of  (1)  before  us, 
let  us  consider  the  case  of  axisyinmetric  waves  on  the  surface  of 
a  shallow  spherical  ocean  of  an  Incompressible,  Inviscid,  gravi- 
tating fluid  of  constant  density.   Let  the  components  of  velocity 
of  a  particle  be  u-,,  the  velocity  along  a  meridian  from  north  to 
south;  Up,  the  radial  velocity;  and  u_  the  velocity  along  a 
circle  of  latitude  from  east  to  west.   Then  the  momentum  equa- 
tions in  spherical  coordinates  (p,9,<j>)  are 

p 

2U-.        U-.    3Ut               3u,             u~,          pu,         UtU^        u-   cot   Q  T 
1    ,   _1  1   ,          1   ,          3 1    .      1  2         3               _        1    gp 

at  pa©  2  9p  p    sin   Q  Wf  p       "  p  "   5p   39 

!^  +  1^  !!^  ,.  „     £^4.        ^3  ^^2        u|       u|  _  ^  ^ 

3t  p     39  2   9p  p    sin   9  d^  P     "     P         "^   "  "5  3p 

3u-        u,    3u,  3U-,  u-  9u-  u^Up        U-u,    cot   P 

9t  p    39  2  ^p  p   sm  9       J^  p  p 


5p   sin  9 


where  p  is  the  pressure,  5  is  the  constant  density  and  g  is  a 
constant  gravitational  acceleration  due  to  a  body  force  directed 
toward  the  center  of  the  sphere.   The  continuity  equation  is 

P29     p    3p       p       psin9  9^ 

Let  us  assume  that  the  motion  does  not  depend  on  the  longitude, 
and  that  the  linear  shallow  water  theory  is  applicable.   This 
theory  is  here  characterized  by  the  assumptions  that  the  radial 
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acceleration  can  be  neglected  and  that  the  non  linear  terms  in 
the  equations  can  be  neglected.   Under  these  assumptions  the 
third  of  the  momentum  equations  Is  automatically  satisfied  and 
the  first  two  reduce  to 


5u,       T 

^^'  at     op  aQ 

and  the  continuity  equation  reduces  to 

? 

T    3Ut  sin  ©   T  3p  Uo 

^^'  sin  ©  9©        +  p  ap      ^* 


In  addition,  we  suppose  that  Up  =  0  at  p  =  a,  the  bottom  of  the 
ocean,;  and  that  p  =  0  at  the  free  surface  whose  equation  Is 

p  -  a  -  h  -  \!/(©,t)  =  0 

where  h  is  the  depth  of  the  ocean  and  \I/(Q,t)  is  the  elevation  of 
the  free  surface  above  the  equilibriiom  position.  The  linearized 
kinematic   condition  at  the  free  surface  is  then 

(7)  u^  =  vl/t(©,t). 

Prom  (5)  we  find 

(8)  p  =  g6(a+h+ij/-p), 
from  which 


to]  ^  =   0.5   ^ 

and  hence   from   (ij.) 
(10) 


1^-   ££i 
at   "  p  ^0  • 


Differentiation  of    (6)   with  respect   to   the   time   gives 


1        ^     3in  o        1  .  i  ^  p2  42  ^  0 


sin  ©  ae  ^"''  -  at        p   ap  ^     at 

which  by  using    (10)    becomes 

ft         d-  a(p^Up. 

-   0.-5   Q  ^  sin  ©  \I/q  +  ^rr ^^  =  0. 

sm   ©  (3©  ^9        5p 

If  v;e  integrate  this  with  respect  to  p  from  p  =  a  to  the  free 
surface  we  find 

"  iin""©^  ^  ^^^  ®  ^Q  +  (a+h+^|/)2u2t(a+h+\l/,e,t)  =  0. 

Then  if  v/e  use  (7)  and  retain  only  first  order  terms  we  see  that 
\1;  must  satisfy 


^11)  -  ifcr©  3©  ^i^  ^  a© -^  ^   ;^2    '     -  °- 


The  energy  equation  associated  with 


c"   3  --•-  e  al  .  2_Ji  =  0 


sm 


sin  ©  a©  ""''  ^   a©   3^2 


■rU 
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can  be  found  by  multiplying  it  by  \\i.    sin  Q   and  integrating  with 
respect  to  0  and  t:- 

71   t  71   t 

°^  ]   J   ^t  #  ^i^  '^  H  ^*^®  -  J  J   Vtt  ^*^^  =  0 


O  "O  CO 


Integration  gives 


-.71  7C 


oM  * 


o 


|(0,t)  sin  0  d©  +  1   \l/^(e,t)  sin  6  d© 


TC  ^7t 


=  c^  \  \!/|(@,0)  sin  e  d©  +  I  \|/|(Q,0)  sin  ©  d©. 
>^  o  o 

The  interpretation  of  this  energy  equation  depends  of  course  on 
the  meaning  of  \{/  and  c.  It  is  most  easily  interpreted  if  i]/  de- 
notes the  radial  displacement  of  a  thin  elastic  spherical  mem- 

2 
brane.  For  this  case  c   is  proportional  to  the  uniform  tension 

and  the  equation  simply  states  that  the  kinetic  energy  plus  the 

potential  energy  of  the  membrane  is  constant. 


V7.J. 


Li'-r 


:i   e.: 
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3«      Bernoulli-type   Solution. 
We  proceed  now   to   solve 

(12)  ^li^&^i^  ®l!=^^tt 


-■jt<0<ix  ;        0<t<(x»      ;        c=  ^^^ 


a 


subject  to  the  following  Initial  conditions.  We  assume  that  i|^ 
is  initially  in  the  equilibrium  position,  that  is 

(13)  »Ke,o)  =  0. 

vie  also  assume  that  the  initial  velocity  in  an  axlally  symmetric 
neighborhood  of  a  point  P,  wMch  we  may  take  to  be  the  north  pole, 
©  =  0,  is  given  by 


Uk) 


0      <      Q      <      0-r       <     % 


e,  <  ©  <  II. 


iiThere,  in  the  neighborhood  -Q      <  ©  <  ©   f(©)  =  f(-©).   In  addi- 
tion, v;e  require  the  solution  to  be  bounded  at  the  poles  vjhen 
f(©)  is  bounded  at  ©  =  0  and  ©  =  -n:. 

If  we  introduce  the  independent  variable 

X   =  cos   © 
the   equation  becordes 


'r^     >      K      ~- 


>«=    (0,P).t 


(C.^:> 


/«-■•' 


Ss    ?>     j.«;-. 


-  r 


Loo 


eXd'i- 


TC    ^    Vi!     ■ 

'.'nei7©5aj: 


-»       —        ^•- 


I     >      K    >      I      9  u,^^       !r':- 
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where 

■;^(x,t)  =  \l/(arc  cos  x,t) 
with  the  initial  conditions 

V(x,0)  =  0 
(16) 

/4.(x,0)  =  f(arc  cos  x)  =  P(x) , 

The  Bernoulli-type  representation  of  the  solution  of  ( l5) , 
(16)  as  an  eigenf unction  expansion  can  easily  be  found  by  using 
the  Legendre  transform 


1 


7L-"'(n,t)  =  I   P^(x)X(x,t)dx 


where  n  is  an  integer  and 


'■n<'^>  =  ^  h  £h'=^'-^>" 


is  the  polynomial  solution  of 


Ird-^^'  dl^n""  +n(n-M)P„=0. 


It  is  well  known  that  these  polynomials  are  orthogonal: - 
1 


Pn(x)P^(x)dx  =  0      n  ?^  ra 


-1 


^n^^)^^     =   J?! 


-1 
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and  that  for  a  piecewise  smooth  arbitrary  function  Xn(x)  pre- 
scribed In  -1  <  X  <  1 

T  ,   00  .,  ' 

y3_(x)  =  i[X],(x+0)  +  X^(x-0)]  =  ^  EI  (2n+l)Xi(n)P^(x)  . 

The  application  of  the  Legendre  transform  reduces  (l5)  and  (16) 
to 

(17)  X'J'^  +  c^n(n+l)%'''  =  0 

(18)  ;i^"""(n,0)  =  0 

.1 


Xt(n,0)  =   j  P(x)P^(x)dx  =  a^. 


-1 

Th©  solution  of  (17),  (l8)  for  n  ^  0  is 

a  sin  ct  yn(n+l) 
X'"*(n,t)  =  ^ 


c  \/n(n+l) 
and  for  n  =  0 


X-'''(0,t)  =  a^t. 


After  inversion  we  have 

^1 

y(x,t)  =  I  Pq(x)  I  P(x)PQ(x)d 

^-1 


(2n+l)P  (x)  )  P(x)P  (x)dx'sin  ct  v/n(n+l) 
,    CO    .  J  T 


^°  n=l  s/n(n+l) 


™»-»   .-V 


^ 


In  terms  of  the  original  variables  this  Is 


.,,._..,  J' 


(19)     (©,t)  =  I  Pq(cos  0)  J   f (9)Pq(cos  ©)  sin  QdQ 

0 

pi 
(2n+l)P    (cos   ©)  I       f (e)P^(cos   Q)   sin  ©dP.sin  ct  /n(n+l) 

1   0°.       ""      -"© 

+  _±_  \  ■ —  II  II      ... 

^°  n=l  yn(n+l) 

The  representation  (19),  although  useful  for  some  purposes, 
does  not  present  a  direct  plcutre  of  the  waves  as  they  travel 
from  the  north  to  the  south  pole.   For  example,  the  wave  front 
requires  some  time  to  get  to  the  south  pole,  and  hence  \J/('n;,t) 
must  be  zero  for  sufficiently  small  t  but  it  is  not  easy  to  see 
that 

e 

t 


^  Pq(-1)  j   f (G)Pq(cos  9)    sin  @d© 


2 


(2n+l)P^(-l)  I   f (9)P^(cos  0)    sin  ©de.sln  ct^  /n(n+lj 


Pr^  C 


n=l  /n(n+l) 

is  equal  to  zero  if  t,  >  0  is  small.   This  remark  of  course 
merely  vinderlines  the  well  known  disadvantage  of  expressing 
solutions  of  wave  motion  problems  as  an  infinite  series  of 
functions  involving  the  time.   We  can  expect  a  clearer  plctvire 
of  the  wave  motion  if  we  can  express  \I/(o,t)  as  the  sum  of  a 
finite  number  of  terms.   The  series  (19)  can  be  siommed  by 
applying  the  Laplace  transform  with  respect  to  t,  and  sub- 
sequently using  various  identities  and  expansions  Involving 
the  Legendre  fixnctions.   Instead  of  doing  this  we  prefer  to 
apply  the  Laplace  directly  to  (15)  and  (16), 


f     ;^  i- 


'r'.4'.n/' 


;?    -^it-''   -u^V 


T"- 


f-'Jt'o  rr'o    ^;''' 


'  ^- f.. 
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U.   D'Alembert-type  Solution. 
Let  us  denote  the  Laplace  transform  of  X(x,t)  by 


oo 


y(x,s)  =  {       e'^^   /^(x,t)dt. 
^o 

The  result  of  applying  this  transform  to  (l5)  and  (16)  Is  easily 
seen  to  be 

and  we  require  the  solution  of  this  equation  to  be  bounded  at 
X  =  -1,  X  =  +1  when  this  is  the  case  for  P(x).   Let  us  set 

so  that  we  have 

V  =  -  1  i  iX 

where 

A  = 


_./:?:? 


2c 

As  we  will  see  it  is  sufficient  to  take 


V  =  -  i  +  U, 


Vie  now  have 


(20)         ^(1-^^)  H  +  v(v+l)X=  -  -^P(x). 


<f)' 


The  homogeneous  equation 


^  (1-^')    g  -^  ^(^■'■l)^  =   ° 


is    satisfied   by   the   linearly   independent   solutions 


16 


Pv(x)    =   P    ^ 


Pv(-^) 


./2 


(x)    =  ^   cosh  An 


n 


^  cosh  An 


,00 

o 

.00 


cos   Azdz 
/cosha+x 


cos   Azdz. 
\/coshz-x 


These  special  Legendre  functions  are  sometimes  referred  to  as 
conical  functions.   The  function  P  (x)  is  analytic  in  -1  <  x  <  1 
but  it  has  a  logarithmic  singularity  at  x  =  -1;  P  (-x)  has  a 
logarithmic  singularity  at  x  =  1. 

In  terms  of  P  (x)  and  ?   (-x.)    it  can  be  shown  that  if  P(x) 
is  piecewise  continuous  in  -1  <  x  <  1  then 


^X 


^(X,S)  =  -^ 


2c  coshAn 


P^(x)  J   P^(-^)F(^)d?  +  P^(-x)|  P^(^)P(^)d^ 
-1  '  X  - 


is  the  only  bounded  solution  of  (20).   In  terms  of  0  the  solu- 
tion is 

.It 


(21)  /(©,s)  =  — ^ 


2c  coshAi: 


P  (cos  ©) \  P  (-cos  a)f(a)sin  a  da 


9 


Q 


_+P  (-cos  e)\  P  (cos  a)f(a)sin  a  da_ 


/  -r  4  T. 


i? 


The  problem  novj  is  to  find  the  Inverse  of  (21)  in  a  form  suit- 
able for  our  purposes.  A  representation  analogous  with  the 
D'Alembert  representation  for  ^^;aves  traveling  along  a  string  can 
be  found  by  performing  the  inversion  after  l/coshA-n  has  been 
expanded  in  powers  of  c"   .  VJe  proceed  to  illustrate  this  by 
investigating  a  particular  case. 

Let  us  consider  in  some  detail  the  case  for  which 


f(©)  =  • 


k  (cons't.)    0  <  ^  <  ©o 


1° 


Q  <  e  <  -IT, 

O       — 


For  this  initial  distribution  of  \I/.  (0,0)  we  have 

re  e 

\  P  (cos  a)f(a)sin  ada  =  k  I  \,i(©  -a)sin  a  P  (cos  a)da 
^0  ^0^ 


where  v-.(^)    Is    the  unit   function  defined   by 

0  <  ^ 

'0  K  <  0 


If  we  use   the   representation 


a 

•D   / « ^ »   ^  \   _  \/A  \        cos  h  Audu 
P„  ( c  o  s    a )    =  ^ \       ; — 

^'-.  /cos   u-cos   a 


valid  for  0  <  a  <  ii,   we   find 


r-v: 


18 


°     e^^sln  u  du 

714/cOS     U-C03     © 


0  >  e 


P  (cos  a)f(a)sin  ada 


_k__ 


-0 


e   sin  u  du 


X/cos  u-cos  9 


0  «  0  <  9   • 
~  o 


-e 


Also 


n 


P  (-cos  a)f(a)sin  ada  =  k 


11 


n(Q^-a)sin  a  P^(-cos  a)da 


6 


and  by  using 


-It 


P  (-cos  a) 


yg  \   cos  h  A(u--iT)du 

^  J  /cos  a  -  cos  u 
a 


we  find 


7t 


0 

-71-©, 


©  >  ©. 


(  P  (-cos  a)f(a)sin  ada 
^© 


k^ 
/In 


e  sin  u  du 


-(1:.©^) 


VCOS  ©+C03  U 


.It-© 


°   e'^^sin  u  du  [ 


y/S-ri    \  VCOS     ©    +OOSU 


'^  <  ^o- 


.(Tt-©^) 


Hence  for   0  <  P  <  0^:- 


rr,    . 


OZi    Ji! 


'v 


•ivt.:-  J.~/^^/ft       '    ■■ 


H     '^  A  .**; 


f*  ''■       '^■'■i' 


»    «   >•  « ) 


»   '■ 


1  ! 


Mi    : 


4  .§--•  1  •' 


-.:  '     O     >     si 
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/P^(cos    P)\ 


e^^sln  u  du     "X 


tj/(©,s)    =  — ^ 


/  ry\     VCOS     ©+COS    U 

-(71  -9; 


2(/2c  cosh  .Alt 


r^-®o      Au  . 
-P  (cos  e)\        e_^i2^^ 


OS  u   \ 


■  9 


+p  (-cos  e) 


Au  .     •, 

e   sin  u  du 


-©    a/cos  u-cos  9 


and  for  9  >  9  :- 
o 


^(©,s)  = 


k 


2/20  cosh  At: 


.P„(-cos  9) 


e   sin  u  du 
-9  A/aos  u-cos  9 


For  simplicity  let  us  confine  our  attention  to  the  nature  of 
the  distxirbance  at  the  north  pole  9  =  0,  which  is  the  center  of 
the  initial  disturbance,  and  at  the  antipodal  point  9  =  ir. 

Since  P_(l)  =  1  we  see  that  the  Laplace  transform  of  the 
disturbance  at  0  =  n  is 


?(7t  ,s)  = 


^9    , 

r  o     Au    .       J 

e   sm  u  du 


2/2c  cosh  Air 


Ai/cos  u-cos  9, 


-9. 


If  we  expand  l/cosh  An  in  poxirers  of  e~   we  have 


k 


00 


(22)     ^(7t,s)  =  -zzIZ  (-1) 

\/2o  m=0 


n 


^o  ^-[(2n+l)7t-u]A_. 


SI 


-9. 


c>Vcos  u  -  cos  9 


n  u  du     ^/ 


I  i  f 


■« ; 
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Now 


^-[(2n+l)7i-u]X 


i:^ 


(2nlt+it-u)    I   ^     0" 

5 J=    --f 


00 


0 


c    /.2      (2n+7i:-u)' 
^     t 


dt 


where  the  Bessel  ftinctlon 


^  q)^"^ 


satisfies 


I    (z)   =  ^ 

^  m=0  raiP{[i+m+l) 


2dy,„dY        /„2,     2v      _^ 


Therefore  the  inversion  of  (22)  gives 

/2c    0 


,,  CO  „ 

(23)    M/(7i,t)  ==  -—n  (-1) 


sin  u 


i._Q     /cos  u  -  cos   © 


•l/[t  - 


(2mt+7t-u) 


] 


•0 


c    \  .2      /2mt+7t-u\ 

iJ^  -  ^     s — ^ 


2 


du  , 


For  any  finite  t  the  series  in  (23)  contains  only  a  finite 
number  of  terms;  and  consequently  it  displays  some  aspects  of 
the  nature  of  the  disturbance  at  ©  =  tc  more  clearly  than  the 


-.0. 


I.      —f,-  .;»■  .-V-.    : 


a<^j';^'^i;irl  £. 
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infinite  series  representations  (19).  For  example  if 
0  <  t  <  — r —  we  can  see  from  (23)  that 

\|/(Tl,t)  =0       0  <  t  <  -rr^ 


so  there  is  no  disturbance  at  the  south  pole  during  this  time 

Interval,  The  front  of  the  disturbance  from  the  circle  of  co- 

latitude  ©  =  ©^  does  not  arrive  at  ©  =  it  until  t  =  ,  During 

the  time  interval  — - —  <  t  <  3  -  - ■  the  disturbance  arrives  at 

c   —   —  -^  c 

©  =  11  from  the  initially  disturbed  area  and  is  then  reflected 
back.   In  this  time  interval  the  disturbance  is  determined  by 

^(^,t)=-^r°   ^in  u        ,,tt,(I^)3l  [^/t^,(SzU)^]du 
l/2c  \    /cos  u-cos  o^  *^     c    O  d^      c 

^'o 


or 


■»o 


y/2C  /cos     U-COS     9^        O    ^  ^  I'         \    c     ' 


M  °  it-9^       It-©, 


°  <  t  <  3^ 


where  M  is  the  maximum  of  —  -  t;  -  ©  ,  By  using  the  last 

°        7t-9° 

integral  it  can  be  shovm  that  if   .    <  t  <  —  the  disturbance  is 
positive  and  increasing  (k  >  0) ;  but  to  see  the  precise  varia- 

71-9 

tion  in      .    °   <  t   <  —  and   the  remainder  of   the    interval 
t        — rt     —  c 

It- 9         ^"0 

<  t  <  3"-r —  ^Q   would  have  to  evaluate  the  integral  numeri- 

~"c      -   c 
cally. 


)•:;  J  o&:,v.iiy 


iJ 
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The  Laplace  transform  of  the  disturbance  at  ©  =  0  is 


r  ri^ 


k 


e      sm  u  du 


?k/l+COSU 


^(0,3)     =  ^       2    ^'^ < 

2/?  c     cosh  X-K    j 


rlt-© 


°        e^^sin  u  du 


or 


(21;)    i|/(0,3)    =  "2 


Ji^ 


L(7:-9  )Vco3   e^+cos   u 
27t-e_ 


k 


\ 


c^d+l+A)^        2/2 c^  cosh  Alt  \ 


e^^^-^^sin  u  du 
Q  Vcos   9q-cos   u 


If  we   expand   1/ cosh  At:  in  poiijers   of  e'  ^   [2k.)   becomes 


00 


27t-9 


?(0,s)    =  iL.  .  ^     ^  (.1) 
s         (/2  c  n-0 


n 


Q, 


°  e-(2ntc4-u)A3.^  u  du 
cA/cos    ©q-COS    u 


The   inverse   of   this   is 


00 


(25)    \l'(0,t)    =  kt   -  ^     ^  (-1) 

v/2c  n=0 


n 


211-9^ 


sin  u 


J  Q  /cOS    Q    -COS    u     V. 


^.[t-(^^iSlii)]. 


/t: 


'■"■0^2    /^2    .2n7t+Uv^ 


e 


-2. 


For  0  <  t   <  -^;    (25)    reduces   to 


>l/(0,t)    =  kt 


6 


;  h    f  f 


'  <-     o'^TC 


•  -     It    ',iy-«- 


\X- 
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which  shows  that  the  disturbance  grows  linearly  with  the  time 

during  the  interval  0  <  t  <  -— ,   In  this  interval  the  distur- 

bance  is  sustained  by  the  arrival  of  disturbances  coming  from 

the  circles  of  latitude  in  the  range  0  <  ©  <  9^,      In  order  to 

discuss  the  change  in  the  disturbance  during  subsequent  finite 

intervals  of  time  we  need  to  evaluate  a  finite  number  of  inte- 

Q                  2ti+© 
grals  given  by  (25).   During  the  interval  — 2  <  t  < for 

example,  the  disturbance  is  determined  by 


sm  u 


\/cos  Q 


h[t.f ]i  [|/t274H^ 

-COS  u  <^    c   O  ^1     ^d 


c  —   ~   c 


It  is  of  some  interest  to  consider  the  disturbance  caused 
by  an  initial  concentration  of  velocity  at  9  =  0.   If  we  take 

K 


(26)        f(©)  =  /  :§         0  <  0  <  e 


>  e  <  ©  <  7t 

we  see  from  (21)  that  the  Laplace  transform  of  the  disturbance 
for  ©  >  e  is 

f^         k 

\j/(©,s)  =  — —^ •  P„(-cos  0)  \  P„(cos  a)  -%   sin  ada. 

2c^cosh  At:    ^        ,   ^        e^ 


:t  lol 


„  rt,;+  .f.  f      '• 


i„] 


';!  \'\ 


fn   zc: 


2k 


The   limit   of   this   as    e  — >  0  is 

(27)  t{e,s)    =  —^ •    P„(-cos   ©). 

1|C  cosh  .Xti 

It  should  be  observed  that  if  (26)  is  substituted  in  the 
energy  equation  of  section  2;  and  e  alloi^/ed  to  approach  zero  we 
find  that  the  energy  associated  with  the  delta  distribution  of 
initial  velocity  is  infinite. 

Since 


7t-9 

P  (.cos  0)  =  iZI   f     ^°^^  ^^ dz     0  <  0  <  71 

V  n      ] 


)q      /^^ 


9   +COS  z 


the  transform  (27)  can  be  written  as 

(ZQ)         ^(0,s)  =  — ^2 1      gp^fcl  ^z   _  ^^ 

l|c  cosh  .Atc   \     \/cos  Q+cos  z 
^  0 


where  A.  =  -•  l/s  -r— .   If  we  express  cos   Az/cosh  Xit  as  an  infinite 
sum  of  partial  fractions  and  then   invert,  we  find  the  Bernoulli- 
type  solution  which  agrees  vjith  (19)  evaluated  for  the  special 
delta  distribution  defined  by  (26)  as  e  — s»0.   The  D'Alembert 
type  of  representation  can  be  found  by  expanding  l/cosh  Xit   into 
a  sum  of  exponentials  prior  to  performing  the  inversion. 
The  equation  (28)  is  equivalent  to 


■ay'     < 


•\ 
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% 


'^ 


71- e 


\     (sin  a)\i;(a,s)da  =  — -^ 
A  2c  cosh  Xot  Jq 


f  / 


cos  ©  +  COS  z  cosh  Xz  dz 


(29) 


Cq/^ 


kV2  e 


-xti  71- e 


2c2(l+e-2^^) 


I 


Xz 


e   /cos  ©  +  cos  z  dz 


-(oi-e) 


=  V!  ^  (_,,n  r"'   ^-[  (2n+l).-z]X  ^,,3  ,  ,  ,,3  ,  ,, 


2c  n=0 


-(7t-Q) 


Since 


-a/s''-b    r  -St 
e        ~  \  ® 

4) 


5(t-a)+  vj(t-a).ab 


I^(bJt^I?)'' 


Jt2. 


2  a2 


dt 


where  5(z)  is  the  Dirac  delta  function,  the  inversion  of  (29) 
gives 


at 


71- @ 


I 


k  j/Z  00 
(sin  a)\l/(a,t)da  =  --^  ^  (-1)" 
e  2c   n=0       ..^^.gj 


^     6(t.[i^2ll)Zz£]j 


+  l^  j^t  -  [i^S±li2Ill]  I .  [i2n+l)7t::zj 


'     T  f  c      .2      r(2n+l)7i-Zi2\  -t 
j  ^H      r  (2n+l)7i-z.,k^ 


/cos   8  +  cos   z   dz 


or 


r\ 


.•*  * 


■yj 
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ct-2mi-© 


9  ^  ct-2(n+lhr+e^ 


(30)    (sin  a)il/(a,t)da  =  ^^  2i:(-l)  ]^ 


,c,.   ^/l[|y^T2tZJ)] 


+  Vj(^)^(t  -^) —  V  •/cos  9  -cos  c(t-^)  d^  . 


For  finite  t  the  series  on  the  right  hand  side  of  (30) 

Q 

contains   only  a  finite  niimber  of  terras.      If  t   <  —  no   terms 
•^  c 

appear  on  the  right,  which  corresponds  to  the  fact  that  the 

Q 

disturbance  does  not  manifest  itself  at  Q  until  t  =  — .   Let  us 

c 

now  consider  in  some  detail  the  time  interval  0  <  Q  <  ct  <  I+ti+O, 
which  covers  the  time  required  for  at  least  two  transits  of  the 
disturbance  at  any  9.  For  this  interval  only  ti-jo  terms  appear 
on  the  right  of  ( 30)  namely 


r%  r   Ct-e  :''ct-27t-9 


ko/2  [  'c  '      k^/2 
(sin  a)il/(a,t)  =  ^    \  -  ^ 


9  J  ct-2:i)9        ^ct-ij.71+9 

c  c 


(31) 

■^1  —  /y  ( 2 1  -^ )  ]  ( 

6(^)  +  n(?)|(t"-?)  ^       '  ■  (  ^cos  9 -cos  c(t  -^;d?. 

^  "^  /?(2t-^)     ) 


For  0<Q<ct<27t-9,  which  is  the  time  interval  during 
which  the  disturbance  pass  ess  from  9  to  the  south  pole  and  back 
to  9,  (31)  yields 


..   '^o 


■  o-  ©■ 


•'.I.' 


x,.-i- 


.,.,    t 


\     (sin  a)\|/(a,t)da  =  —^r-  /cos   e-cos   ct 


ct-Q 


+     Ol_  (t-^)  -iJ ycos   e  -cos   ctt-^)    dK 


from  which 

(32)        *(e,t)  = 


k 


(/2 


licy'cos   Q-cos   ct 

k^/H  r*  ^         i^(fPr?) 


8  \ 

e/c 


8         \  ycos   Q  -  cos   cZ       -f  ,  2  _     2 


Since  \lf(9,t)  =  0  for  ct  <  0,  (32)  shows  that  as  t  passes  through 
the  value  t  =  ©/c  the  disturbance  at  ©  suffers  an  infinite  jump 
discontinuity.   If  ct  is  slightly  larger  than  ©,  ij;(©,t)  behaves 
essentially  like  const,/  yet-©.   For  0  <  ©  <  ct  <  271  -  ©,  (32) 
gives  the  wake  of  the  disturbance  left  behind  at  ©  while  the 
front  travels  to  the  south  pole  and  is  there  reflected  back. 

After  the  second  passage  of  the  front  through  ©  while  the 
front  passes  to  the  initial  center  of  the  disturbance  and  then 
back  to  the  position  ©,  that  is  for  211  -  ©  <  ct  <  Su  +  ©,  the 
disturbance  at  ©  is  given  by 


:<& 


'T 


CT:!' 


r-T  . 


11 


sin  a  ^(a,t)da 


28 


ct-e 


ico^a  \    = 


ct-27t+e 


or 


2Tc-e 


(33)    *(e,t)  =  -g—  \ 

O/c 


/cos  ©-COS  cZ^t^-i^^ 


d?  . 


A  comparison  of  this  with  (32)  shows  that  the  disturbance  ex- 
hibits the  same  kind  of  singular  behavior  at  ct  =  271-©  as  it 
does  at  ct  =  ©. 

After  the  third  passage  of  the  wave  over  the  position  © 
while  it  passes  on  to  the  antipodal  point  to  be  again  reflected 
back,  that  is  for  27t  +  Q  <  ct  <  i;?:  -  e  the  "tail"  or  wake  of  the 
disturbance  at  6  is  given  by 


It 

6 


sin  a  \Ka,t)da 


k.^/2 


-ct-9 
c 


ct-2Ti+e 


^^"^^  I,  [S/^(2t-^T]  yc(3s  e  -cos  c(t-^)  d? 


oT /cos  0-COS  ct 

2c   " 

,    _  /Ct-27I-© 

ko/2 


—^^--^—  i.i^y/zUt^K)]  /COS  e-cos  c(t-^;d^ 


•A- 


\r 
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or 

kQ/2 
(3ij-)    ^(<?,t)   =  -  — 

[(.C/COS    ©-COS    ct 


°        \  /cos   e-cos   c^  /l2~32 

^0/c  ^  '^ 


kQ,/2    ^ 


^id  /?^  d^ 


\/cos    Q-cos   c^  /,2   yZ 
27t+Q  '^ 


c 


This  shows  that  the  intensity  of  the  singular  part  of  the  dis- 
turbance at  ct  =  271+6  is  the  same  as  it  was  at  ct  =  211-©  and 
ct  =  a  although  there  is  a  reversal  of  sign.   The  general  ex- 
pression (30)  shows  that  the  singular  part  of  the  front  of  the 
disturbance  changes  only  in  sign  and  these  changes  occur  at 
passages  3,5,7  •••   •   The  maintenance  of  the  strength  of  the 
singularity  is  not  remarkable  because  under  the  assumend  initial 
conditions  the  disturbance  is  charged  ^^:ith  infinite  energy, 

VJe  see  now  that  lire  can  proceed  to  find,  and  evaluate 
n\imerically  if  need  be,  the  various  integrals  which  give  the 
disturbance  at  any  position  0  <  Q  <  n  for  any  interval  of  time. 

The  disturbance  at  the  north  pole  during  the  interval 
0  <  ct  <  1;tc  can  be  found  by  allovjing  9  to  approach  zero  in  the 
integrals  (32),  (33)  and  (3i|).   The  wake  left  at  ©  =  0  as  the 


^-rM 


^^/ 


4v    <D. 
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disturbance  emanates  from  the  north  pole,  proceeds  to  the  south 
pole  and  comes  back  (0  <  ct  <  Zn)    is,  as  we  can  see  from  (32), 


(35)  >l'(0,t)  = 


k 


0 


ct: 
tpc  sin  -g- 


IT  I  sin  c^   — :z=ir  a^' 


0 


sm  c^ 
2 


v^^^^:? 


Notice  here  that  if  ct  approaches  Sir  from  below  >|/(0,t)  behaves 
like 


*(°»*^  ^      2^^  +  b  ^(2at-ct). 


Equation  (3^1-)  can  be  written 


(36)    ^(e,t)    =  - 


ko/^ 


l].c  ycos  ©  -  cos   ct 


\J 


-2-11-9 


©/c 


a.  (§•/?!?) 


■1^2 


sin  ^  yriTin^  I  csc^  f/^' 


dK 


^0  r 


V 


^^(l  \/?^) 


2it+© 


an  f /l-sin2  I  csc2  f /^I? 


dK   . 


The  change  in  sign  of  the  second  integral  is  due  to  the  fact 
that  /cos  Q-cos  c^  is  positive  and  hence 


/cose-cos  c^  =  -/Z   sin  ^  /1-sin  ^  esc  ^  when  27t  <  c^  <  iiu. 

The  limit  of  (36)  as  ©  — >  0,  which  gives  the  disturbance  at 
9  =  0  for  the  interval  2it  <  ct  <  I^ii,  is 
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2      l)  ^   (sin  T5-J  •*  t  -^ 


where  the  principal  value  of  the  integral  is  to  be  taken.  If 
ct  approaches  Zii  from  above,  v;e  see  from  (37)  that  \l/(0,t)  be- 
haves like 


ij/(0,t)  ^      _  ^^^  +  b^(ct-2Ti). 


Zir-ct 


The  last  results  show  that  as  the  front  of  the  disturbance 
arrives  at  the  initial  center  of  disturbance  at  e  =  0  it  ex- 
hibits a  more  violent  behavior  than  it  does  when  it  is  away 
from  the  north  pole  at  the  position  ©  where  0  <  ©  <  tc« 

The  disturbance  at  the  south  pole  for  the  interval 
0  <  ct  <  l].7i  can  be  found  by  taking  the  limits  of  (32),  (33)  and 
{2h-)    as  ©  approaches  it.  Equivalent  results,  hox-jever,  can  be 
found  with  less  trouble  by  substituting  ©  =  71:  directly  in  (27). 
If  we  do  this  we  find  [since  P  (1)  =  1]  that 

Ttk 
^(it,s)  =      " 


^c  cosh  '^'^ 


(38) 


k^oi  00 


=  10|  ^.  (.i)n^-(2n-.l),rX 
2c^  "T" 


fto' 


ri-^o    ;•-;' 
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(We  cannot  find  \|r(0,t)  by  directly  substituting  9  =  0  in  (2?) 
because  P  (-cos  Q)   has  a  logarithmic  singularity  at  ©  =  0,   The 
inversion  of  (38)  gives 


(39)   ^(Jt,t)  = 


Tlk-   03 


5[t  -  (2n+l)5] 


2c  n=0 


(-1)''  J 


"N 


+  ^^  [t-  (2n+l)f].(2n+l)|  '\   . 


1/ 


f"  -  (2n+l) 


2  n' 


c 


L    yt2-(2n+l)2^ 


2^ 


c 


This  shows  that  if  0  <  ct  <  7i  then  ^in,t)   =   0,      If  ti:  <  ct  <  3ii 
we  have 


\l'(7t,t)  = 


irk 


2c 


|8(t.f)+|v|(t-f) 


^  ^i(  ai*  "r?  ) 


/ 


t   — ^ 

c 


The  delta  function  singularity  which  the  front  exhibits  as  it 
arrives  at  the  south  pole  at  ct  =  t:  is  stronger  than  what  it 
has  as  it  arrives  at  ©  (0  <  ©  <  ii)  but  loss  strong  than  the 
singularity  exhibited  when  the  front  comes  to  the  north  pole. 
As  the  disturbance  recedes  frora  the  south  pole 
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This  Is  the  wake  left  at  0  =  Tt  diorlng  the  Interval  r.  <  ct  <  311 
while  the  front  of  the  distiirbance  moves  from  south  pole  to 
north  pole  to  south  pole  again.  For  the  interval  3it  <  ct  <  5^^ 
we  have 
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In  this  interval  the  disturbance,  as  it  comes  back  to  9  =  ii, 
exhibits  except  for  sign  the  same  kind  of  singularity  it 
possessed  at  its  first  arrival.   The  wake  left  dviring  the 
interval  371  <  ct  <  5*^  Is 
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